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We study SU(2) gluodynamics at finite temperature on both sides of the deconfining phase tran- 
sition. We create the lattice ensembles using the tree-level tadpole-improved Symanzik action. 
The Neuberger overlap Dirac operator is used to determine the following three aspects of vac- 
uum structure: (i) The topological susceptibility is evaluated at various temperatures across the 
phase transition, (it) the overlap fermion spectral density is determined and found to depend on 
the Polyakov loop above the phase transition and (Hi) the corresponding localization properties of 
low-lying eigenmodes are investigated. Finally, we compare with zero temperature results. 
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I. INTRODUCTION 

More than ten years ago, in a model generalizing random matrix theory, M. A. Stephanov [1[ has 
built in the feature that in St/ (3) gluodynamics above T c the different Polyakov loop sectors should 
behave differently, fn the sectors with complex-valued average Polyakov loop the chiral condensate 
was expected to turn to zero at a temperature T substantially higher than T C) the deconfinement 
temperature. Prior to this paper, similar observations in quenched 517(3) lattice simulations had 
been reported by Chandrasekharan and Christ [2|, Q- For SU{2) lattice gluodynamics, where the 
Polyakov loop is real, it was predicted that the chiral condensate stays non-zero, < iptp >^ 0, for all 
temperatures T > T c in the sector with a negative average Polyakov loop, L < 0. 

Simulating quenched SU(3), Gattringer et al. [J] came to a different result. These authors used the 
spectral gap in the spectrum of the Dirac operator as an order parameter for the restoration of chiral 
symmetry. It was found that the spectral gap opens up at one single temperature T — T c in all three 
Z 3 sectors. 

In SU(2) gluodynamics this question has not been clarified before. In this case the behavior may 
be completely different. Indeed, there are pure vacuum gauge field configurations with L = 1 for 
which the Dirac operator has trivial zero modes for periodic fermionic boundary conditions. Contrary 
to SU(3), for SU(2) (as well as for any even number of colors) through a non-periodic Z(2) gauge 
transformation with —Is Z(2) this is equivalent to the existence of zero modes for antiperiodic fermion 
boundary conditions in the case of L = — I . In the light of this argument, we shall examine whether 
Stephanov's prediction for the Dirac spectrum remains valid in the case of SU(2) gluodynamics in the 
deconfined phase. Preliminary results were presented in [5|. 



We base our study on Neuberger's overlap Dirac operator [6j]. Our interest will be concentrated on 
its spectrum, how it is related to the average Polyakov loop, and how this is reflected by the localization 
properties of the modes in various parts of the spectrum. Because of the index theorem [7], |8(, the zero 
modes give us easy access to the topological charge Q of the lattice configurations, which is sufficient 
to calculate the topological susceptibility. In our present work, in contrast to Ref. [9(, we will not 
exploit the capability of the overlap Dirac operator to provide us with a tool to analyse the topological 
charge density [10| . 

A first study of the overlap Dirac operator D ov in SU(2) gluodynamics at finite temperature has 



been undertaken by the authors of Ref. [ll|. They have used the usual Wilson gauge field action 
which means that for the same lattice spacing their configurations are substantially more rough than 
ours. We compare our results with theirs at several places in the following. 

The outline of the paper is as follows. In Sect. [IT] we explain the action used and quote the ensembles 
investigated. In Sect. Mil we discuss the structure and handling of Neuberger's overlap operator. In 
Sect. [TV] we report about the T-dependence of the topological susceptibility and in Sect. |V] about 
the temperature and Polyakov loop dependence of the spectral density. In Sect. |VTJ we define the 
spectral gap as a deconfincment order parameter. In Sect. IVHI the localization properties of the 
Dirac eigenfunctions are described in the confined and deconfincd phase, and the dependence on the 
Polyakov loop is pointed out for the latter. In Sect. IVIIII we assess the volume and lattice-spacing 
dependence of all discussed observables. In a separate Sect. lIXI we consider an exceptional, relatively 
well separated part of the non-zero modes with very small eigenvalues. After a discussion of topological 
susceptibility, spectral density and localization for the case of zero temperature in Sect. |Xl a summary 
is given in Sect. IXI1 



II. THE IMPROVED ACTION 

In an initial study, partly reported at Lattice 2007 Q, we have analyzed ensembles of 0(100) 
statistically independent SU{2) configurations, generated in the quenched approximation with the 
tadpole improved Symanzik action on lattices of size 20 3 x 6. This action is known to suppress 
dislocations, that would lead to unphysical near-to-zero modes of the Wilson-Dirac operator. These 
are highly unwanted because they lead to difficulties to choose a uniform p parameter (see below). 
The form of the action is : 

S = Amp 2^ Spl - ^| 2_, S rt ' W 

pi ° rt 

where S p i and S r t denote the plaquette and 1x2 rectangular loop terms in the action, S p ^ r t = 
(1/2) Tr (1 — U p i : rt). The factor u = (Wixi) 1 ^ 4 is the input tadpole factor. It is selfconsistently 
determined from W± x i = ((1/2) Tr U p i) computed at zero temperature [12]. The deconfining phase 
transition, which is of second order for SU(2) gluodynamics, is not easy to locate precisely [l3| because 
of strong finite volume effects. We have approximately determined the critical /3i mp , c = 3.248(2) on 
lattices of size 20 3 x 6, i.e. for N T — 6. This corresponds to a ratio T c /y/a = 0.71(2) [9J. 

To check for finite spatial volume effects, we made additional simulations also on 24 3 x 6 lattices as 
well as on smaller lattices 16 3 x 6 and 12 3 x 6. Furthermore, we also checked whether our results are 
close to the continuum limit. For this check we have simulated the neighbourhood of the transition on 
a finer lattice, with N T = 8, for which a critical /3i mp . c = 3.425(5) was found in simulations on lattices 
with a spatial extent N s = 24. For these 24 3 x 8 lattices our actual measurements were performed 
only at one temperature very close to T c in the confinement phase. 

The parameters of the lattices used in this study and the respective statistics are reported in Table [Ij 
The lattice scale a needed to compute T/T c has been obtained by interpolating the results of Ref. [12| • 

For a comparison with zero temperature, we have made simulations on symmetric lattices. Re- 
spective parameters are also given in Tabic |TJ At zero temperature the inverse coupling was chosen 
according to the lattice size such that the physical size was kept approximately fixed and equal to 
(1.4 fm) 4 . 



N T x iV s 3 


Pimp 


T/T c 


#conf. 


■* * modes 


6x 12 3 
6x 12 3 


3.20 
3.23 


0.91 
0.97 


998, 910 
893, 900 


20 
20 


6x 16 3 
6x 16 3 
6x 16 3 
6x 16 3 


3.20 
3.23 
3.325 
3.50 


0.91 
0.97 
1.15 
1.55 


482, 517 
246, 262 
183, 215 
308, 96 


30 
50 
30 
50 


6x 20 3 
6x 20 3 
6x 20 3 
6x 20 3 
6x 20 3 
6x 20 3 
6x 20 3 
6x 20 3 
6x 20 3 


3.20 
3.23 
3.275 
3.30 
3.325 
3.35 
3.40 
3.50 
3.64 


0.91 
0.97 
1.05 
1.10 
1.15 
1.20 
1.31 
1.55 
2.00 


48, 49 
152, 147 
89, 102 
42, 56 
101, 98 
123, 75 
148, 
89, 98 
99, 98 


50 
50 
50 
50 
50 
50 
50 
50 
50 


6x 24 3 
6x 24 3 
6x 24 3 
6x 24 3 


3.23 
3.325 
3.35 
3.50 


0.97 
1.15 
1.20 
1.55 


55, 41 
115, 50 
46, 00 

78, 81 


45 
45 
45 
45 


8 x 24 3 


3.42 


0.99 


28, 31 


25 


10 x 10 3 
12 x 12 3 
14 x 14 3 
16 x 16 3 
20 x 20 3 


3.096 
3.196 
3.281 
3.3555 
3.5 









200 
100 
100 
100 
99 


40 
40 
40 
40 
20 



TABLE I: The simulations details. The fourth column shows numbers of configurations with L > and L < 
separately where the two states can be distinguished. 



III. THE MASSLESS OVERLAP DIRAC OPERATOR 



The massless overlap Dirac operator has the form 6] 



D ov ='-[ 1 
a 



^ 



(2) 



where D\y = M — p/a is the Wilson-Dirac operator with a negative mass term p/a, M is the Wilson 
hopping term with r — 1, and a is the lattice spacing. With respect to the locality property of the 
overlap operator the optimal value of the p parameter is found to be p — 1.4 also for the lattice 
ensembles under investigation. Standard boundary conditions, anti-periodic in time and periodic in 
spatial directions, are imposed to the fermionic field and the overlap Dirac operator. 
In order to compute the sign function in 



D 



w i 



D^Dw 



75 sign(H w ) 



(3) 



where Hw = 75 Dw is the hermitian Wilson-Dirac operator, we have used the minmax polynomial 
approximation. More precisely, we have treated 20 to 50 lowest Wilson-Dirac eigenmodes explicitely 
and used the polynomial approximation in the subspace orthogonal to these lowest Wilson-Dirac 
eigenmodes. The overlap Dirac operator D ov constructed this way preserves the chiral symmetry 



even on a lattice with finite a and allows to exploit the perfect chiral properties of the emerging Dirac 
overlap eigenmodcs. The overlap operator will be called D in the following and replaces the continuum 
Dirac operator D = D^ 7^, where D^ — d^ — igA^ is the covariant partial derivative with the gauge 
field background A^. 

IV. THE TOPOLOGICAL SUSCEPTIBILITY xto P (T) 

We solved the Dirac equation numerically finding its eigensystem 

D ip„ — A„ ip n , (4) 

restricted to the 50 lowest eigenvectors by means of ARPACK [14| . They are perfectly located on the 
Ginsparg- Wilson circle. By stereographic projection they are mapped to the imaginary axis. These 
purely imaginary eigenvalues Ai mpr are then representing the 0{a) improved overlap Dirac operator. 
The spectral densities presented later are distributions with respect to the imaginary part of Ai mpr . 

Our first task was the search for exact zero modes. Their number is related to the total topological 
charge Qtop of the lattice configuration through the index theorem [3, Q : 

Qtop = Qindex = N- - N + , (5) 

where N- and N+ are the numbers of zero modes ifiio with negative and positive chirality, tp i0 75 ipio = 
±1, respectively. Actually, we found only configurations with N + or A_ or both vanishing. For the 
lattice ensembles the average topological charge (Qtop) should vanish, while the second moment {Ql op ) 
measures the strength of global topological fluctuations. The topological susceptibility is measured as 



(Q 



2 
top/ 



Xtop = —^— , (6) 

where V is the four-dimensional lattice volume in physical units. In Fig. [T] (left) we show a histogram 
of the topological charge in the confinement phase, close to the transition. Fig. Q] (right) shows the 
corresponding histogram for a temperature higher up in the deconfincment phase. The histogram in 
the confinement phase can be approximately fitted by a Gaussian distribution, in the deconfinement 
phase the Gaussian fit is of more poor quality. Let us now discuss the topological susceptibility as 
function of temperature. In Fig. \2\ (left) we show the topological susceptibility xtop separately for 
the positive (L > 0) and the negative Polyakov loop sector (L < 0) over the full temperature range. 
At temperatures T < 1.31 T c we have obtained both sectors in the course of the usual Monte Carlo 
sampling. Flips between the sectors happened even at T > T c due to finiteness of our volume. For 
T = 1.55 T c and 2.0 T c we have not observed such flips and stayed within the L > sector. We then 
generated configurations for the L < sector in the following way. We applied a flip transformation 
changing the sign of the Polyakov loop to an equilibrium configuration with L > and then used the 
Monte Carlo algorithm to produce a necessary number of configurations for the other sector. No flips 
back to the L > sector have been observed during these runs. Taking all temperatures into account, 
there is no systematic influence of the sign of the Polyakov loop or, what is equivalent, of the type 
of temporal boundary conditions, on the topological susceptibility. This conforms to the observation 
in Ref. [15j . for quenched SU(3) configurations at T = selected to have topological charge Q = ±1, 
that the content of zero modes does not change with a continuous, complex- valued factor used to 
implement non-trivial temporal boundary conditions in the Dirac operator [311 ] - 

The measurements for the two signs of the Polyakov loop agree with each other at all T within two 
standard deviations. In Fig. [5] (right) we compare our final averages for Xto P (T), which include for 
T <Tp all configurations and for T > T c only the subsample with L > 0, with the results of Alles et 
al. [Ty] [32|]. These authors have presented the values of 10~ 4 x Xtopl 'A* for various values of (3 for 
Wilson's action representing different temperatures. We took A^ = 14.15(42) MeV [Ty] and extracted 
their susceptibility values Xtop(T) from these data. The topological susceptibility is slowly decreasing 
with increasing temperature for both sets of data. Notice that the overlap definition of Q results in a 
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FIG. 1: Probability distributions of the topological charge Q for two temperatures below and above T c on the 
20 3 x 6 lattice. 
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FIG. 2: The fourth root of the topological susceptibility xto P as function of T on the 20 3 x 6 lattice. Left: 
separately for L > and L < 0. Right: comparison of our final result with that of Alles et al. [JU| and Edwards 
et al. fill. 



systematically higher susceptibility than the improved field theoretic definition employed by the Pisa 
group. Additionally we show results obtained in Ref. [Il(. We took the data for the topological charge 
from their Table IV for temperatures T/T c = 1.0, 1.4, 2.0 and computed the respective susceptibilities 
Statistical errors were computed by the bootstrap method. One can see that the results from Ref. 
agree well with ours for T/T c = 1 and T/T c = 1.4 but slightly disagree for T/T c = 2.0. 
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FIG. 3: The third root (ivp(X)/V) 1 ' 3 of the spectral density of non-zero eigenmodes of the overlap Dirac 
operator (at A = approaching the third root of the quark condensate) is shown for two temperatures T < T c 
as measured on the 20 3 x 6 lattice. The modes are counted with a bin size of 30 MeV. 



THE SPECTRAL DENSITY AND CHIRAL SYMMETRY RESTORATION IN 

DIFFERENT Z 2 SECTORS 



The chiral condensate (ipip) is related to the density p(X) of the non-zero eigenvalues A at A 
via the Banks-Casher [17| relation: 



w> 



lim lim lim 

A^0 rn^O V-nx 



V 



(7) 



The non-zero modes are globally non-chiral, but still locally chiral over a wide part of the spectrum, and 
their pseudoscalar density is known to be correlated with lumps of the topological charge density [18| . 
In the chirally broken phase the required limit |0 of p(X) is non- vanishing at A = [l7|. In the 
chirally symmetric phase we expected vanishing p(X) = in a finite region around the origin, in other 
words, that the spectrum develops a gap. It cannot be excluded that the gap (an operative definition 
is given below) tends to zero in the infinite-volume limit while the spectral density at zero remains 
vanishing, p(X = 0) = 0. For the confinement (chirally broken) phase we find indeed that the spectral 
density in physical units is practically T independent and slightly increasing towards A = 0. This can 
be seen in Fig. [3] where the third root (np(X)/V) 1 ^ 3 of the spectral density is shown. At A = it 
should approach the third root of the quark condensate. 
The eigenvalue density is defined as 



P(A) 



lim 

AA^O 



(N(X,AX)) 
AA 



(8) 



where N(X, AA) counts the number of eigenvalues per configuration falling into the bin intervall 
[A — AA/2, A + AA/2]. The averaging over all configurations available is denoted by (...). In practice 
we have chosen a finite bin size AA that is mentioned in each case. 

We shall note that in the confinement phase, comparing results for configurations with average 
Polyakov loop L > and L < 0, we found at low A the quantity (ttp(X)/V) 1/3 to be 50 ~ 70 MeV 
higher for the negative Polyakov loop sector. We observed that this difference comes mostly from 
configurations with large \L\. This observation leads us to believe that such difference disappears in 
the thermodynamic limit since it is known that |_L| — s- in this limit. 

For the deconfinement phase, when we take only configurations with an average Polyakov loop 
L > 0, the left panel of Fig. Q] shows that p(X) non-uniformly decreases with increasing temperature, 
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FIG. 4: The third root (7rp(A)/V) ' 3 of the spectral density of non-zero eigenmodes of the overlap Dirac 
operator shown for 4 temperatures T > T c as measured on a 20 3 x 6 lattice, evaluated separately according 
to the sign of the averaged Polyakov loop: L > (left) and L < (right). The modes are counted for L > 
with a bin size of 50 MeV, for L < with a bin size of 30 MeV. 



indicating the decrease of the quark condensate until a gap in the spectrum opens at T/T c = 1.55 and 
gets wider at T/T c = 2. This does not completely agree with our expectation above. The gap, when it 
finally opens (for the L > sector) at T/T c rs 1.5 needs a careful investigation in the infinite- volume 
limit. For configurations with L < the right panel of Fig. [4] shows that the third root of np(X)/V at 
low A stays non-zero and even grows with increasing temperature. We checked that decreasing of the 
bin size down to 10 MeV does not uncover any gaps in the spectrum for L < 0. 

For the temperatures above T c in the sector with L > 0, in the interval A < 50 MeV we observe a 
particular enhancement in the spectrum of non-zero modes, that monotonously decreases with rising 
temperature. We will refer to these modes as "near-zero modes" . A similar observation has been 
made by Edwards et al. 11] for quenched configurations generated with the Wilson action. The fact 
that it reappears with the improved gauge action strongly hints that it is not a lattice artefact. The 
ensemble average of the lowest-lying one of these non-zero eigenvalues seems to be a natural definition 
of the spectral gap (see Sect. IVip . The presence of the spectral enhancement, however, suggests a 
special, more careful analysis with respect to its origin that will be given in Sect. IIXI and will result 
in an alternative definition of the spectral gap. 

The sector with L > plays a particular role because, as our results indicate in agreement with 
model considerations of Ref. [1], in this sector, together with the standard antiperiodic boundary 
conditions for fermions, chiral symmetry can be restored at sufficiently high temperature. We will 
furthermore see on the opposite that for the sector L < 0, with antiperiodic boundary conditions kept, 
chiral symmetry is not restored in the sense of the Banks-Cashcr relation. 

However, with dynamical fermions included the sector with L > is the dynamically chosen sector 
in the high-temperature phase. For quenched configurations, when the fermionic boundary conditions 
can be changed at will from antiperiodic to periodic boundary conditions in time, the role of the 
sectors L > and L < is simultaneously interchanged. 



VI. THE SPECTRAL GAP 



The spectral gap g\ has been defined as the ensemble average of the smallest non-zero eigenvalue. 
In Ref. [4[ Gattringer et al. have shown for SU(3) gluodynamics that the gap, as a function of 
temperature, has a different magnitude but a similar behavior for the real and both complex sectors 
(defined according to the complex phase of the averaged Polyakov loop) . The phase transition - defined 
by the respective g\ becoming non-zero - occurs at the same T c . With increasing lattice volume the 
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FIG. 5: The spectral gap for SU(2) lattice gauge theory as function of temperature, evaluated separately 
according to the sign of the averaged Polyakov loop on a 20 3 x 6 lattice, L > (left) and L < (right). Notice 
the scales for the gap differing by a factor 10. 



gap g\ at each temperature T > T c has the tendency to decrease, but this decrease seemed to stop 
at a finite limit. The infinite- volume limit of the gap in the three sectors has not yet been carefully 
analysed. 

SU(2) gluodynamics has only two sectors in the deconfincmcnt phase, distinguished by the sign 
of the (real-valued) averaged Polyakov loop. We show in Fig. [5] (left) a clearly defined and rapidly 
growing gap g\ for configurations with L > 0, whereas for configurations with L < the gap remains 
very small up to temperatures several times higher than T c , as can be seen from the right panel of 

Fig.m 

Our data presented in section ITXl (see Fig. flTT right)) strongly suggest that the small gap for L < 
is a pure finite- volume effect and vanishes in the limit of spatial V3 — > 00. 

Physically more important in the light of our expectations would be to prove that the gap for L > 
- although also decreasing with increasing volume - has a finite limit at infinite volume. We postpone 
the discussion of this question to Sect. IIXI Although the separation of the near-zero modes from 
the definition of the gap, that we will discuss there, will make the existence of a finite gap for all 
temperatures above T c and for all available volumes obvious, a real infinite-volume limit of the gap 
g\ has still to be explored. 



VII. LOCALIZATION IN DIFFERENT PARTS OF THE SPECTRUM 



The scalar density of an eigenmode i[)\(x) corresponding to an eigenvalue A is denoted as s\(x) — 
ip x (x)ip\(x), such that ^2 s\(x) — 1 by virtue of normalization. The inverse participation ratio (IPR) 



I\ is the natural measure of the localization. For any finite volume V it is defined by 



h=Vj2sl 



(x). 



0) 



The IPR characterizes the inverse volume fraction of sites forming the support of s\ (x) — ipl {x)ip\ ( x ) ■ 
Its value is equal to the volume V for modes which are localized on a site, scales with V if it is localized 
on a finite number of lattice sites and is independent of the volume if its localization volume grows 
with the lattice volume. For T = the localization of overlap eigenmodes has been first investigated 
in Refs. pj, [H, HI H3 , for T ^ in Ref. [H, SI- With a chirally improved Dirac operator, the 
localization properties have been investi gate d in detail earlier by Gattringer et al., for quenched 
SU(3) Yang-Mills theory at T = in Ref. [24[ and for finite temperature in Ref. [25|. The localization 
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FIG. 6: The IPR averaged over zero modes and over subsequent spectral bins of width 30 MeV for two 
temperatures T < T c on the lattice 20 3 x 6. 



L>Q 



L<0 



yo 

80 
70 
60 
50 
40 
30 
20 
10 



I ' T/T c = 1 .05 — 
"f 1.20 — — 1 
1.55 — ^ - 

2.00 — 

* 

• • * * * * 




12 

10 

8 

> 6 

4 
2 



T/T c = 1.05 ^+- 

I 1.20 — — 

I 1.55 — ^- 

. ] 2.00 — : 


200 400 600 800 1000 1200 14 
X, MeV 


00 


100 200 300 400 5( 
X, MeV 



FIG. 7: The IPR averaged over zero modes and over subsequent spectral bins of width 70 MeV (left) and 
30 MeV (right), separated according to the sign of the averaged Polyakov loop L > (left) and L < (right). 
Results are shown for four temperatures T > T c on the lattice 20 x 6. 



properties and their change with temperature are surprisingly similar for SU(3) and SU(2). This 
applies to the overlap operator as well as to the chirally improved Dirac operator. 

In the following figures the IPR are shown at A = averaged over the zero modes, followed by 
the IPR averaged over the subsequent bins of non-zero modes. From Fig. [6] we conclude that for the 
temperature near but below T c the IPR (localization) increases with decreasing eigenvalue, more or 
less in a monotonous fashion. There is no clear mobility edge. The zero modes are even more localized 
exceeding the first bin by w 50 %. Thus, out of the low-lying modes, the higher ones are continuously 
less localized. We found that at these temperatures for configurations with large in absolute value 
negative Polyakov loop, L < 0, the modes are less localized by a factor 2^3. Again, as for p(A), this 
difference should disappear in the thermodynamic limit. 

In Fig. [7] (left) we show for L > that with increasing temperature in the deconfinement phase 
the average IPR within the respective eigenvalue bins is increasing. With increasing temperature, the 
effect moves to higher and higher eigenvalues A. This can be considered as a mobility edge moving, 
together with the gap, to larger A in the deconfinement phase for L > 0. In the negative Polyakov 
loop sector the IPR is taking values at a much lower level. The temperature dependence at fixed A 



10 

is reverse compared to the L > sector, and the temperature dependence becomes stronger towards 
small A. The spreading of the IPR values with respect to temperature is almost twice as large for zero 
modes compared to the smallest non-zero modes. 



VIII. FINITE VOLUME AND FINITE LATTICE SPACING EFFECTS ON 
TOPOLOGICAL SUSCEPTIBILITY, SPECTRAL DENSITY AND LOCALIZATION 



240 



220 



© 200 



180 



160 



140 




0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 
T/T„ 



FIG. 8: The fourth root of the topological susceptibility Xto P as function of T on 12 3 x 6, 16 3 x 6, 20 3 x 6, 
24 3 x 6 and on 24 3 x 8 lattices compared with the results of Ref. [la ]. 
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FIG. 9: The volume dependence of the third root of the averaged spectral density of non-zero eigenmodes 
of the overlap Dirac operator at T = 0.97 T c (left) and T = 1.15 T c (right). While N T — 6 in all cases, we 
compare four spatial lattice sizes with N s — 12, 16, 20 and 24 below T c and three spatial lattice sizes with 
N s — 16, 20 and 24 above T c . The modes are counted with a bin size of 30 MeV below T c and with a bin size 
of 50 MeV above T c . 



To study finite volume effects we have additionally generated ensembles of configurations on 12 3 x 6, 
16 3 x 6 and 24 3 x 6 lattices (see Table H]). In Fig. [5] we compare our results for the topological 
susceptibility Xtop obtained on all our lattices. One can see that the finite volume effects are rather 
small and unimportant for this observable. Finite lattice spacing effects seem to be also negligible, 
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FIG. 10: The third root of the spectral density of non-zero eigenmodes of the overlap Dirac operator for two 
temperatures T <T C closely approaching T c on a 20 3 x 6 and the 24 3 x 8 lattice, respectively. The modes are 
counted with a bin size of 30 MeV. 

since the \top result on 24 3 x 8 lattice is in accord with that on the 20 3 x 6 lattice at the temperature 
where we have data available (see Fig. |5J). 

From Fig. [5] (left for T < T c and right for T > T c , restricted to the L > sector) we conclude, 
that the spectral density (if divided by the 4-volume) does not depend on the lattice volume for 
large A, while some volume dependence is seen for A < 100 MeV for T < T c (Fig. [9] left). Fig. [10] 
demonstrates a weak, if any, dependence of the spectral density in the confinement phase when the 
transition temperature is approached. 
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FIG. 11: The IPR averaged over zero modes and over spectral bins of width 30 MeV (left) and 50 MeV (right). 
Left: comparison for T = 0.97 T c (N T = 6) between four spatial volumes with N 3 — 12, 16, 20 and 24. Right: 
comparison for T = 1.15 T c in the L > sector between three spatial volumes, on 16 3 x 6, 20 3 x 6 and 24 3 x 6 
lattices. 



Coming now to the localization properties of the modes, we conclude from Fig. [TT] that the IPR 
with increasing spatial volume increases (comparing lattice sizes from 12 3 x 6 to 24 3 x 6). This applies 
both to the confinement (left) and the deconfmed (right) phases. In Fig. 1 121 we show, comparing the 
20 3 x 6 with the 24 3 x 8 lattice (at approximately the same temperature and with a ratio of physical 
three-dimensional volumes of 0.86) that the lattice spacing effects are small. However, this is not the 
case for the zero modes. Their IPR nearly doubles when the lattice spacing a decreases by a factor 
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FIG. 12: The IPR averaged over zero modes and over spectral bins ol width 30 MeV checked for lattice spacing 
effects by comparison between the 20 3 x 6 lattice and the 24 3 x 8 lattice for two temperatures T < T c close to 
each other. 



6/8. This corresponds to a dimensionality of the zero modes close to d sa 2 [2C 



IX. COMPARISON WITH EDWARDS ET AL. 



As it has been said before, the overlap operator spectrum at finite temperature has been studied 
both in SU(2) and ST/ (3) gluodynamics by Edwards et al. [ll|. In SU(2) case they made simulations 
at T/T c = 1.0, 1.4 and 2.0 . It has been found that for all temperatures considered the spectrum 
(for antiperiodic boundary conditions and positive Polyakov loop) consisted of two parts: a group 
of unexpected (and unobserved with staggered fermions !) near-zero modes, i.e. modes with small 
eigenvalues below 0.05/a with a spectral density decreasing as a function of A, and the rest with larger 
eigenvalues above 0.05/a with a spectral density increa sing with increasing A. There was a large gap 
found between these two parts. The authors of Ref. [ll[ assumed that the small eigenvalues (the 
near-zero modes) are due to instanton-antiinstanton pairs and presented some evidence supporting 
this assumption. We prefer to call them more generally "topological objects" in what follows. In 
particular, Edwards et al. found that the number of small (near-zero) modes n sm , combined with 
the number of zero modes no, is equal to the number of level crossings in the spectral flow of Hw, 
pointing towards a common origin. They stated that (Q 2 ) agreed well with the average (n) where 
n = n sm + no- We will find a similar remarkable agreement. Furthermore, it was found that the 
distribution of n — n sm + jiq is in good agreement with a Poisson distribution 



e -(n) 
P(n;(rc)) = — («)" 

77' 



(10) 



where (n) is the ensemble average. 

Our Fig. Q] (left) demonstrates that we have also observed the two parts of the eigenvalue spectrum 
mentioned above, although these two parts are not separated by a gap for T/T c < 1.55, rather by a 
plateau of low spectral density. We therefore defined a boundary between the two parts of the spectrum 
for T/T c < 1.55 at 150 MeV which is approximately located at the minimum of the distribution p(A). 
We checked that all our results discussed below vary only within error bars when we vary this border 
value by 50 MeV. 

Following the discussion in Ref. [Il| we studied properties of (n) and of the probability distribution 
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FIG. 13: The density of topological objects (see text) in units of T* as a function of temperature T, from our 
data obtained for two 3- volumes at T > T c and from the data of Ref. [Til l. 



P(n; (n)) defined above. In Fig.[T3]we present the density p top defined as 



Ptop ~ y 



(11) 



as a function of the temperature. For comparison, we show data from Ref. [1 II ] which we took from 
their Table II. We computed the respective statistical errors by applying the bootstrap method to 
data collected there. We present in Fig. Q2] our results for the 24 3 x 6 lattice as well. There is good 
agreement between our results and results of Ref. [11| for smaller temperatures while data clearly 
disagree for the largest temperature T/T c = 2.0 [33|. One can see no essential volume dependence for 
this quantity which has been concluded by Edwards et al. as well. It is also seen that the density pt op 
is slowly decreasing with increasing temperature. The decrease with temperature is similar to that 
of the topological susceptibility. Numerically we indeed found nice agreement between x arid ptop as 
can be seen from Fig. [T3J 
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FIG. 14: Comparison of p\' /T c (full symbols) and x 1/4 /T c (empty symbols) for 20 3 x 6 (squares) and 24 3 x i 
(circles) lattices. 
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In Fig. [T5]the ratio (n)/D is depicted, where D is the variance of the multiplicity n of topological 
objects including near-zero and zero modes. For the Poisson distribution exactly (n) = D holds. One 
can see that our results agree with this equality except for one data point. Thus we confirm the 
earlier observation of Ref. [Tl[ that the distribution of n is consistent with the assumption of mutual 
independence between the respective topological objects. 

In Fig.[16]we show the ratio (n sm )/(n) which indicates how the contribution of the near-zero (small 
eigenvalue) modes to all topological objects changes with volume and eventually with temperature. 
At two temperatures, T/T c — 1.15 and T/T c = 1.55, where we are in the position to compare two 
three-dimensional volumes (N s — 20 and N s — 24 with the same N T = 6) we could find this ratio 
rising with the volume. In other words, the larger the three-dimensional volume, the more near-zero 
modes exist compared to the zero modes. This behavior is consistent with a binomial distribution 
proposed in Ref. 11] for topological objects with positive and negative charge, with no = \n+ — ?i-| 
zero modes and n sm = 2 min(n+,n— ) near-zero modes. On the other hand, the share of near-zero 
modes among all topological objects seems to decrease with increasing temperature. One should have 
in mind that the physical spatial volume in our setting rapidly decreases (like (N S /(N T T)) ) with 
increasing temperature. For T/T c > 1.2 the ratio is already less than one half, indicating that the 
volumes are already so small that the total multiplicity is suppressed such that there are more zero 
modes than near-zero modes. 



Q 

V 




FIG. 15: The ratio of the average number of all topological objects (n) to its variance as a function of 
temperature T, from our data obtained for two 3- volumes at T > T c and from the data of Ref. [ill ]. 



Now let us look at the problem of the spectral gap taking into account the existence of a separate 
part of near-zero eigenvalues in the spectrum as discussed above. Alternatively to the definition of 
the spectral gap used in section IVll one can consider a definition g\ which takes into account only the 
bulk non-zero modes in the spectrum, thus ignoring the small non-zero eigenvalues (near-zero modes) 
altogether. This alternative definition can be applied to our data at T/T c = 1.55 , a temperature at 
which we observe a clear gap between near-zero and larger eigenvalues. In the left panel of Fig. [T7] 
we show the dependence of the conventional spectral gap, as defined in section I VII on the inverse 
lattice size for two temperatures in the case of L > 0. The general tendency of a spectral gap is 
decreasing with increasing volume. The interesting question is whether it has an infinite- volume limit 
different from zero. One can see that for T/T c — 1.15 the spectral gap decreases very steeply with 
increasing volume, starting from values like g\ « 300...400MeV on lattices of size L ~ 1.4fm. The 
gap is already less than 100 MeV for L > 2.2 fm, eventually turning to zero in the infinite-volume 
limit. At higher temperature, T/T c — 1.55, we also see that the "old" gap g\ decreases with increasing 
volume, although not so fast. Still, we cannot exclude that the fate of the spectral gap is the same at 
the higher temperature. 



Thus, the old definition of the spectral gap gives some support to the conclusion made in Ref. [11 1 
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FIG. 16: The ratio of the average number of small non-zero modes (n sm ) to the average number of all 
topological objects (n) as a function of temperature T, from our data obtained for two 3-volumes at T > T c 
and from the data of Ref. \Vm . 
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FIG. 17: The spectral gap g\ as defined in section [Vl] vs. inverse lattice size for T/T c = 1.15 and T/T c = 1.55, 
separated according to the sign of the averaged Polyakov loop L > (left) and L < (right). 



that in the deconfinement phase of SU(2) gluodynamics the chiral condensate is nonzero, at least up 
to some specific temperature, even in the L > sector. 

In the right panel of Fig.[T7]the conventional spectral gap g\ is also shown for L < 0. Compared to 
the case of L > 0, the spectral gap for L < is much smaller. There is a clear tendency for g\ at L < 
steeply to decrease with increasing volume. Data in Fig. Upright) suggest that the dependence of g\ 
on 1/L is 1/L 1 with 7 > 1. Our rough estimation results in 7 ss 3. To make a more precise statement 
we need new data points for bigger volumes. The same is valid for the data set at T/T c = 1.15 for 
L > shown in the right panel of Fig. Q~7] 

For the same physical 3-volume, the gap increases with increasing temperature only in the case of 
L > 0. For the negative Polyakov loop sector the tendency is reverse, corresponding to the increasing 
spectral density. 

Let us now turn to Fig. [18] where results for the "new" definition g\ of the spectral gap, excluding 
the small eigenvalues of the near-zero modes, are presented. Now we see a strong indication for a 
convergence of the spectral gap g\ to a finite intercept at 1/L = for both temperatures. Note that 
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FIG. 18: The "new" spectral gap g\ for L > 0, as defined in section HXl with exclusion of the near-zero modes, 
is shown vs. inverse lattice size for T/T c = 1.15 and T/T c = 1.55. 

for T/T c = 1.15 only eigenvalues A < 120 MeV have been excluded in this way, a value substantially 
below the value of the new emerging gap. Thus, excluding the near-zero modes leads to a definition 
of the spectral gap g\, such that it stays non-zero in the limit V — > oo for L > at both temperatures 
T>T C . 

Summarizing our findings presented in this section we confirm the existence, at T > T c , of particular 
modes with small non-zero eigenvalues which are most probably related to topological objects, similar 
to the zero modes. So it is natural to consider them separately from the other (bulk) non-zero modes. 
We are planning to enter a further study of their properties in the near future. 



X. COMPARISON WITH T=0 



The details of simulations at zero temperature are presented in Table HI The lattice coupling and 
the lattice size were chosen in such a way that the physical box size was kept approximately fixed and 
equal to 1.4 fm. In Fig. [111] we show results for the topological susceptibility as a function of the lattice 
spacing squared. The linear fit results in x 1 — 202(5) MeV. This number is in good agreement with 
earlier results obtained with various gluonic definitions of the topological charge on lattices generated 
with the Wilson action. It also agrees well with our results for T <T C . Comparing Fig. fTS] with the 
results for the Wilson action [19( , obtained using the same method to compute the topological charge 
as applied in the present paper, we see that the scaling properties are better for the tadpole-improved 
Symanzik action. 

In Fig.[5nithe spectral density as obtained on the 14 4 lattice for /3i mp = 3.281 is shown. This coupling 
value is close to one value (/3; mp = 3.275) used in the finite temperature simulations at T ~ T c . One 
can see that the spectral density in Fig. [20] is in good agreement with that shown in Fig. [3] except for 
the lowest interval A < 100 MeV. For T = 0, we do not see in this range of A the increase of p(\) with 
A — > that was observed approaching T — > T c from below. Comparing with the respective result for 
the Wilson action [19[ we see good agreement. 

We have also measured the IPR for zero temperature. The results are presented in Fig. [5TJ The 
results are again in full agreement with those of Ref. [19} ■ Comparing the T = results with our 
observations at T < T c (see Fig. [5]) at approximately equal /3i mp , we see an essential difference: at 
non-zero temperature the IPR close to A = is substantially higher. 
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FIG. 19: The fourth root of the topological susceptibility at T — as function of the lattice spacing a. The 
empty circle shows the result of the linear extrapolation to a = 0. 
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FIG. 20: The spectral density obtained on the 14 lattice with j3 — 3.281. The modes are counted with a bin 
size of 50 MeV. 



XL SUMMARY 



After a first overlap study of topological aspects of SU{2) pure gauge theory at finite tempera- 



ture [ll|, based on the Wilson gauge field action, we returned to these topics. We have performed 
measurements of the topological susceptibility, the spectral density and the localization properties of 
fcrmionic eigenmodes with the help of the overlap Dirac operator in finite temperature SU(2) glu- 
odynamics, simulated with the tadpole-improved Symanzik action. We found that the topological 
susceptibility in the confinement phase is almost independent of temperature. It starts slowly to de- 
crease in the deconfinement phase, in agreement with previous results |16| which were obtained with 
the Wilson gauge field action and an improved gluonic definition of the topological charge density. In 
contrast to the other quantities, we did not find systematic effects of the sign of the averaged Polyakov 
loop on the topological susceptibility. 

While these results (see Sect. IIV[) could have been expected, the remaining observations concerning 
the spectral properties might come as a surprise, if not at a qualitative then at a quantitative level. 

Assuming antiperiodic boundary conditions for the fermionic fields throughout this paper we have 
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FIG. 21: The IPR averaged over zero modes and over subsequent spectral bins of width 30 MeV at T = for 
four different /3-values . 



discussed the behavior of the chiral condensate defined via the Banks-Casher relation below and above 
the phase transition. A qualitatively different behavior of the spectral density has been observed at 
T > T c for ensembles having opposite sign of the Polyakov loop L (see Sect. El . The spectral density 
is decreasing with rising temperature for L > until it is finally disrupting around T ~ 1.5 T c . On 
the other hand, for L < the density begins to rise with increasing temperature. 

Although the spectral gap g\, defined as the average of the lowest non-zero eigenvalue, as a function 
of temperature shows qualitative agreement with respective data for 5(7(3) gauge theory (see Fig. [SJ 
we observed that the respective volume dependence is completely different. We found that g\ tends 
fast to zero with increasing volume for L > at T/T c = 1.15 and for L < for both T/T c = 1.15 and 
1.55 in contrast to the situation in SU(3) where a rather weak volume dependence of the gap was 
observed in both real and complex L sectors. 

Increasing of the gap g\ with increasing temperature for L > and fixed lattice size is explained 
by the strong volume dependence of the number of near-zero modes as can be seen from Fig. 1161 
We are making measurements on the lattice, with N s = 20 defining its physical size, which decreases 
with increasing temperature, resulting in a decreasing number of near-zero modes. The consequence 
is that not each lattice configuration has these modes. This in turn makes the ensemble average g\ to 
increase with temperature. 

The localization properties of the eigenmodes have been investigated in terms of the inverse partic- 
ipation ratio (IPR), and the dependence of the average IPR on volume and discretization has been 
monitored all over the lower part of the spectrum (see Sect.EHl We repeat here only two remarkable 
facts. At first, sizeable lattice spacing effects have been found only for the zero modes (close to T c ) 
which is pointing to a dimensionality sw 2. Second, at L < 0, the eigenmodes responsible for the 
non- vanishing chiral condensate above T c , that are filling the region up to A ~ 400 MeV, are very 
delocalized compared to the isolated modes that exist in that spectral range for L > 0. 

Special considerations have been devoted in Sect. IIXI to the near-zero modes found in the sub- 
ensemble with L > 0. The spectral density has been divided into an interval with small eigenvalues and 
the rest with all larger eigenvalues, motivated by the observation that for large enough temperatures 
T rs 1.5 T c these two parts of the spectrum become well separated, whereas a small number of near- 
zero modes still exists. The topological susceptibility (which had been defined through the total 
topological charge Q obtained by counting the zero modes alone) now has been compared with the 
density of topological objects defined including the near-zero modes. Good agreement was found. This 
justifies the dilute-gas approximation, x = Ptop with charge Q rs ±1 objects. This is corroborated 
by the observation that the multiplicity of all topological objects follows the Poisson distribution 
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with reasonable accuracy as indicated by Fig. [T51 and was reported earlier in Ref. [ll[. Further, the 
percentage of near-zero modes among the lower part of the spectrum (i.e. including also the true 
zero-modes) has been measured for different temperatures and volumes, see Fig. [TH] It was found to 
rise with increasing volume, and the apparent temperature suppression observed in our measurements 
might be mainly explained as a volume effect. 

Based on all these observations, emphasizing the special role of the near-zero modes, an alternative 
definition of the spectral gap g\ excluding the near-zero modes has been proposed. Following this 
definition, a non-vanishing infinite-volume limit of g\ in the deconfined phase at all temperatures 
T > T c now seems possible for the L > sub-ensemble, in clear distinction to the other definition 
and, as we have seen, to the L < sub-ensemble. Then we end up with an unusual situation in the 
temperature interval T c < T < 2 T c : a finite density of the near-zero modes coexists with a gap 
between these modes and the bulk of the spectrum. Whether we really may speak about broken chiral 
symmetry depends on a nonvanishing infinite- volume limit of the spectral density p sm (A) at A = 0. 
Otherwise we could ignore the presence of the near-zero modes for this question and consider the new 
definition g\ of the spectral gap as a proper one. In any case, more studies of the thermodynamical 
limit are necessary in order to confirm or reject the possibility to reconcile, in that limit, the restoration 
of chiral symmetry, p(A = 0) = 0, with a gap g\ — ► 0. [3J] 

It shall be noted that by now, the "kinematic" understanding for the relation between chiral sym- 
metry breaking (and restoration) and the sign of an eventually nonzero value of the Polyakov loop (or, 
equivalently, on the change of the boundary conditions in the Dirac operator) is much more developed 
than at the time when our study was begun J5|. 

A series of papers [26|, [27], [28|, [23, (301 discusses the connection between the spectrum of a (fairly 
general) Dirac operator and the breaking of center symmetry i.e. the emergence of a nonvanishing 
Polyakov loop. The most concise quantity in this respect is the "dual quark condensate" defined in 
Ref. [29j , the Fourier transform of the standard quark condensate 

v J = I( Tr [(m + ^)- 1 ]) 

with respect to the angle <j) modifying the boundary condition of the Dirac operator D<p. A "fat 
Polyakov loop" formed out of singly wrapped pathes is obtained as the lowest (n — 1) Fourier compo- 
nent. Higher Fourier components correspond to multiply wrapped Polyakov loops. The condition for 
a nonvanishing "fat" Polyakov loop is a nonvanishing periodic component oc cos(</)) in the fermionic 
spectral density, in other words, a nontrivial response of the spectral density to the change of the 
boundary conditions. 

A microscopic explanation of the onset of such a sensitivity of the Dirac spectrum to the type of 
boundary conditions in the context of a nonvanishing Polyakov loop of either sign needs to be worked 
out. What comes to mind is the interplay of holonomy and topology of dyons [9j. This difference is 
accompanied by a different localization behavior of the lowest fermionic eigenmodes in the two sectors. 
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